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ABSTRACT 



We present Spitzer Space Telescope data of four isolated white dwarfs that were pre- 
viously known to harbor circumstellar gaseous disks. IRAC photometry shows a significant 
infrared excess in all of the systems, SDSS0738+1835, SDSS0845+2257, SDSS1043+0855 
and SDSS1617+1620, indicative of a dusty extension to those disks. The 4.5/im excesses 
seen in SDSS0738, SDSS0845, and SDSS1617 are 7.5, 5.7 and 4.5 times the white dwarf con- 
tribution, respectively. In contrast, in SDSS1043, the measured flux density at 4.5 microns 
is only 1.7 times the white dwarf contribution. 

We compare the measured IR excesses in the systems to models of geometrically thin, 
optically thick disks, and find that we are able to match the measured SEDs to within 3a 
of the uncertainties, although disks with unfeasibly hot inner dust temperatures generally 
provide a better fit than those below the dust sublimation temperature. Possible explanations 
for the dearth of dust around SDSS 1043+0855 are briefly discussed. Including our previous 
study of SDSS1228+1040, all five white dwarfs with gaseous debris disks have significant 
amounts of dust around them. It is evident that gas and dust can coexist around these 
relatively warm, relatively young white dwarfs. 



Introduction 



Since the discovery of the first dusty disk around the white dwarf G29-39 ( Zuckerman &: Becklin 
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to the fate of planetary systems when their host stars undergo post-mai n sequence evolution . 



The most likely origin for the disks is the tidal distruption of asteroids (IGraham et al.lll990 



Debes fc Sigurdssonll2002l : |Jurall2003l ). and the subsequent accretion of this tidally distrupted 
material can explain the large metal abundances in the white dwarf atmospheres. 



The discovery of a gaseous disk around the hot white dwarf SDSS J1228 -I-1040 (IGansicke et al 



2006) and our subsequent identification of a dusty extension to that disk ( iBrinkworth et al. 
20091 ) showed th at debris disks ar e not restricted to low-temperature white dwarfs, as pre- 
viously believed (IKilic et al.l 120061 ) . and dust and gas re adily coexist in suc h a disk. Indeed, 
measured parameters for the 18 dusty white dwarfs (see lFarihi et al.ll2010l . for an overview) 
show that the presence of dust does not appear to correlate with temperature. 

Here we present infrared photometry and SED modelling for a further four white dwarfs 
similar to SDSS J1228+1040: SDSS J104341.53+085558.2, SDSS J084539.17+225728.0, 
SDSS J161717.04+162022.4 and SDSS J073842.57+183509.6 (hereafter SDSS1043, SDSS0845, 
SDSS1617 and SDSS0738, respectively). As with SDSS J1228+1040, all have previously been 
identified from Sloan Digital Sky Survey spectroscopy, reveal i ng Ca II 860 nm emission line s 
from a gaseous disk component (IGansicke et al.l 120071 . |2008| . |2012| ). (IGansicke et al.ll201ll ). 
leading us to search for an infrared exces s indicative of dust y extensions to these disks. 
SDSS0738 was in dependently d i scover ed by lDufour et al.l (I2OIOI ) because of extremely strong 
metal pollution. iDufour et al.l ( I2OIOI ) also reported an infrared excess in the near-IR, but 
did not notice the Call emission of the gaseous disk component. 

SDSS1228, SDSS1043 and SDSS1617 are all hydrogen-dominated DA white dwarfs, 
while SDSS0845 and SDSS0738 are helium-dominated DB white dwarfs. All are metal- 
polluted, all are moderately "warm", with temperatures ranging from 13,000 - 22,000 K, and 
hence they are relatively young (cooling ages of ~100-600 Myr). 



2. Observations and Data Reduction 



2.1. Spitzer Space Telescope 



We were awarded 2.7 hours in Cycle 4 (P40048), 0.7 hours in Cycle 5 (P50118) and 
2.5 hours in Cycle 7 (P70012) to search for a dusty component to the disks around our 
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four targets (see Table[I]for a full list of obse rvations). For SDSS1043, we obtained Infrared 
Array Camera (IRAC; see iFazio et al.l l2005l ) data in all four channels (3.6 /im - 8.0/xm). 
For SDSS0845, we obtained IRAC data at 4.5 /im and 8.0 fim, Infrared Spectrograph (IRS; 
Houck et al. l200i) Blue Peak- Up Imaging at 16 fim and Multiband Imaging Photometer 
(MIPS; iRieke et al.l booi ) imaging at 24 /im. For SDSS1617 and SDSS0738, we obtained 
Warm IRAC data at 3.6 /im and 4.5 /im. 

The IRAC data reduction for all sources was carried out on the individual corrected Basic 
Calibrated Data (CBCD; cbcd.fits) frames downloaded from the Spitzer archive. The CBCDs 
were initia lly combined with the standard SSC software, MOsaic and Point Source Extraction 
(MOPEX: iMakovoz et al.ll2006l ) with dual outlier rejection to create a single mosaicked image 
for each channel. This was used to determine the position of the target. All of the original 
downloaded CBCD frames were then corrected for array location dependence, using the 
correction frames supplied by the Spitzer Science Center (SSC). Aperture photometry was 
carried out on the corrected CBCD frames with IRAF, using an aperture radius of 3.0 pixels. 
The sky subtraction was applied using an annulus from 12-20 pixels in radius. 

SDSS1043 was found to have a blended neighbour in the extraction aperture. We 
calculated the individual flux densities for each of the blended sources by measuring the 
combined flux, then splitting this by the ratio of peak brightness of a Gaussian fit to each 
source in each channel. 

The IRAC photometry for all sources was then converted from MJy sr~^ to /iJy, and 
aperture-corrected using the standard aperture corrections provided by the SSC. We did 
not apply a pixel phase correcti on to the channel 1 data, but included an extra 1% in our 
error budget to account for this (IReach et al.ll2005l ). No colour correction was applied, since 
we quote the isophotal wavelengths, thus eliminating most of the colour dependency of the 
flux calibration and rendering the magnitude of the correction negligible compared to our 
uncertainties. Uncertainties were estimated using the RMS scatter in the photometry in the 
individual frames, divided by the square root of the number of frames. The quoted uncer- 
tainties are e i ther o ur calculated uncertainty or the IRAC calibration uncertainty derived by 
Reach et al.l (120051 ) . whichever was the larger. 



The IRS Peak-Up Imaging and MIPS data for SDSS0845 were combined in much the 
same way, using MOPEX to create a mosaic to measure the source positions, and carrying out 
the photometry on the individual Basic Calibrated Data (BCD) frames. The sky background 
was found to be variable across the IRS Peak-Up array, and so we created a sky flat by 
combining all of the dithers and rejecting the sources. This sky flat was subtracted from 
the individual frames before carrying out the photometry, bringing the sky background to 
zero. The photometry for both the IRS and MIPS data was obtained via IRAF, using an 
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aperture radius of 3.0 pixels with no sky subtraction for IRS, and an aperture radius of 3" 
(1.224 pixels) with a background aperture of 20" - 32" for MIPS. After converting from MJy 
sr~^ to /iJy, the flux densities were aperture-corrected to the calibration apertures using the 
standard corrections provided by the SSC 

Finally, for all of the data for the four targets, we rejected the photometry from the orig- 
inal frames that had been flagged by the MOPEX Dual Outlier rejection algorithm. We take 
the unweighted mean of the remaining photometry in each channel as the quoted flux den- 
sities. An estimate of the uncertainties was measured using (RMS scatter) / sqrt (Number 
of frames), consistently giving uncertainties lower than the instrument calibration uncer- 
tainties. We therefore adopt the instrument calibration uncertainties from the SSC as the 
uncertainties on our data. For MIPS this is 4%, and for IRSPUI this is 5%. 



2.2. AAT, SOFI, UKIDSS, GALEX, SDSS, 2MASS and WISE Data 



Near-infrared imaging of SDSS0845 was obtained at the Australian Astronomical Ob- 
servatory's Anglo- Australian Telescope (AAT) using the InfraRed Imaging Spectrograph 2 
(IRIS2) on 22 February, 2008, and retrieved from the public AAT archive. A 9-point dither 
pattern was used, with individual exposure times of 30s, 42s and 67s, and total exposure 
times of 270s, 378s and 603s in J, H and Kg, respectively. SDSS1617 was observed in 
the Kg band on 1 May, 2010 using SOFI on the ESO New Technology Telescope (NTT). 
A total of 20 5s images were obtained, offset randomly within a 15" box. The data were 
sky-median subtracted, astrometrically calibrated using the WCSToolfl and co-added using 
SWarJ^ Object d etection and the aperture photometry were done using the SExtractor 
(IBertin fc Arnoutd 119961 ). Instrumental magnitudes were converted to the 2MASS system 
using a set of five bright stars near SDSS0845 and ten stars near SDSS1617. The uncertainty 
in the calibration is 0.02 mag in all bands, and was added in quadrature to the statistical 
uncertainties in the photometry. The J, H and Kg magnitudes, fiuxes and uncertainties are 
reported in Table[2J 

Observations of S PSS 1043 in F , J a nd K bands have been obtained with UKIRT/WFCAM 
as part o f UKIDSS Jlawrence et all l2007h on 22 Mav and 15 Dec, 2009 (TableBI). The 



GAL EX (IMartin et al.ll2005f ). SDSS flAbazaiian et all l2009f ) and 2MASS flSkrutskie et al. 



20061 ) archives and the WISE Preliminary Data Release were searched for additional pho- 



^ ht t p : / / tdc-www . harvard . edu/ wcstools 
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tometry of our targets. The archival and WISE photometry are also listed in Table[2j 



2.3. Possible contamination of SDSS J0845+2257 



The NIRI imaging of SDSS0845 obtained by lMelis et al.l (120101 ) revealed the presence of 
a faint galaxy at a 1.8" separation north-east of the white dwarf. We have downloaded the 
NIRI J, if, and Kg-hand images of SDSS0845 from the Gemini Science Archive, and, after 
correction for non-linearity using the python scripts provided by Gemini, processed them 
in the same manner as the AAT and NTT data described above. In order to measure the 
relative contributions of the star and the galaxy, first an annulus was defined, centred on the 
star, with inner and outer radii of 4 and 7" on the sky (large enough to avoid the galaxy) 
to measure the sky background. The star's contribution was then measured from the flux 
summed over a 1.6" diameter aperture, while the galaxy's flux was measured from the union 
of pixels covered by three overlapping circular apertures covering its major axis. To estimate 
the degree to which the galaxy flux was contaminated by light from the star, three identical 
control apertures were placed diametrically opposite to the galaxy apertures relative to the 
star. The instrumental magnitudes were photometrically calibrated using a bright 2MASS 
star in the field. We find J = 16.44 ± 0.05, H = 16.30 ± 0.06 and K = 16.04 ± 0.07 for 
SDSS0845, and J = 20.33 ± 0.79, H = 19.29 ± 0.32 and K = 18.65 ± . 24 fo r the nearby 



galaxy. These results agree well with the NIRI magnitudes of iMelis et al.l (120101 ). apart from 
the fact that we have somewhat larger uncertainties. 

To estimate the contamination of the galaxy in the Spitzer beam, we obtained the 
galaxy spectral energy distributions (SED) of Maraston (1998, 2005) based on the initial 
mass function of Kroupa (2001). We then redshifted all galaxy SEDs by redshifts < z < 3, 
in steps of 0.2, and discarded all combinations of redshift and age that would disagree with 
the observed J — H and H — Kg colours and the apparent Kg magnitude. The NIRI data 
constrains the redshift of the galaxy to 1.4 < -2 < 3.0, with the range 2.4 — 3.0 only feasible 
for young (< 1 Gyr) galaxies. The spatial extension of the galaxy is suggestive of a redshift 
at the lower end of the range allowed by the NIRI photometry (Ferguson et al 2004). 

The contribution of the galaxy to the IRAC-1 and 2 channels is 0.04 ± 0.01 mJy and 
0.05 ± 0.01 mJy, relatively independent of the galaxy type. The maximum contribution in 
the IRAC-3 and 4 channel is ~ 8mJy for a young (< 1 Gyr) high-redshift {z ~ 3) galaxy. 
However, given the combined constraint from morphology and near-infrared photometry, the 
contribution of the galaxy to the IRAC-3 and 4 channels is likely to be < 5mJy. Galaxy 
models are poorly constrained in the mid-IR, and so we are unable to learn anything more 
abotu the galaxy from the IRS-PU and MIPS channels, other than the fact that the galaxy 



- 6 - 



flux continues to drop out ot longer wavelengths. At these longer wavelengths, the SED is 
dominated by the flux density from the dust disk. In Figure[2]we show an 8 Gyr old, redshift 
z=1.6 galaxy as an example of the flux contribution to the Spitzer bands. 



3. Results and Modeling 



We plotted the SEDs of the four WDs, comparing them to models of isolated white 
dwarfs and other system components as described below. We deflne an infrared excess as 
signiflcant if any of the infrared photometric points lies more than 3cr above the model of 
the white dwarf. 

Where we detect an infrared excess, we attempt to model it using an optically t hick, 
geometrically thin dusty disk, wit h a temperature p r ofile T dink oc where /3 = 3/4 (jjura 
20031 ). as previously described in iBrinkworth et al.l ( l2009l ). The ratio of stellar radius to 
distance sets the absolute scale of the white dwarf photospheric flux, and is flxed in our 
modeling based on the best available parameters. The free parameters in the model are 
the inner disk radius (Rm), the outer disk radius {Rout), and the disk inclination (i). A 
grid of disk models was calculated with inner (Tj„) and outer {Tout) disk temperatures, each 
corresponding to a value of and Rout ranging from 100 to 1800 K in steps of 50 K, and 
inclination ranging from to 90 degrees in steps of 5 degrees. Tout "was flxed to be cooler 
than Tin- A least m ethod was used to flt the H—,K— and IRAC-band infrared fluxes. 



See iGirven et al.l (120121 ) for further information on the fltting technique. 



The model is plotted as a dash-dot grey line in the flgures. For each of our white dwarf 
targets, we include two flgures: one showing the best-flt model to the datapoints, regardless 
of physical feasibility, and one showing the closest possible physical model to that best 
flt. Analyses of the photospheric abundances of white dwarfs hosting dusty disks suggest 
that the bulk composition of the circumstell ar material is dorainated by Mg, O, Si, a nd Fe, 
similar to the bulk composition of the Earth ( iZuckerman et al.ll2007l : iKlein et al.ll201l[). The 



condensation temperature of the main minerals carrying these elements is 1400K (ILodders 
20031 ). and we therefore do not expect that dusty debris discs around white dwarfs can survive 
temperatures much above that. Hence, we deflne an inner dust disk temperature of 1400K 
as a physical limit. For each target, we also describe the coolest inner disk temperature that 
we were able to match to the data while staying consistent to within 3a of all the datapoints. 
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3.1. SDSS J0738+1835 



The flux densities from GALEX and Spitzer are given in Table[2l and are plotted in 
light grey in FigurelH alo ng with the SPSS p R7 spectrum (grey line). We also included 
near-IR photometry from iDufour et al.l (l2010l ) , also plotted in dark grey. The DB white 
dwarf model is shown as a dashed black lin e, with Teff = 13600 K, log(g) = 8.5, Rwd = 
0.00886 R© (Table 3; see lDufour et al.ll2010[) . We use a distance of d = 120 pc, rather than 



the 136pc±22 quoted in lDufour et al.l (120101 ). as this gives a better flt to the WD SED. This 
is still consistent with Dufour et al.'s value. 

The photometry of SDSS0738 shows a clear flux density excess in the infrared wave- 
lengths, compared to the white dwarf model. The best fltting disk model (Figured^) has 
Tin = 1800 K, Tout = 830 K and an inclination of 60 degrees, which is consistent with the 
modeling results from iDufour et al.l ( 120101 ). The best physically-motivated model with Tj„ 
= 1400 K, has Tout = 930 K and a face-on inclination of degrees (FiguredJo). This corre- 
sponds to inner and outer disk radii of Rin = 12.4 Rwd (0.11 i?©) and Rout = 21.3 Rwd 
(0.19 Rq). The coolest possible inner disk temperature is constrained by the need for a high 
enough flux density to match the JHK datapoints, and has Tin = 1350 K, Tout = 880 K and 
an incl ination of degrees . These models are both consistent with the modelling carried 
out bv bufour et~aD J2012h . who used Tin = 1600 ±100 K, Tout = 900 +100/-200K and an 
inclination of 40±20 degrees. 



Analysis of the Call emission lines gives an outer radius of th e gaseous disc compo - 
nent of Rout — l-Resin^i, i.e. a maximum value of Rout ~ ^WRwd, (iGansicke et al.ll2012l ). 
While the Call emission lines provide no formal lower limit on Rout, a low-inclination, prac- 
tically face-on geometry seems unlikely, as it would imply a very narrow ring for the gaseous 
component, narrower than the dusty component. 



3.2. SDSS J0845+2257 



The GALEX, AAT and Spitzer photometry for SDSS0845 are given in Tabled and 
are plotted in light gr ey in FigureEl alon g wi th the SPSS pR7 s pectrum, SDSS photometry, 
NIR photometry from lMelis et al.l (120101 ) and lFarihi et al.l (120101 ). and the WISE Preliminary 
Data Release (dark grey). The data show a clear excess at wavelengths longer than 2/im, 
above the flux density expected from the white dwarf and the contaminating background 
galaxy. For the white dwarf contribution, we adopted a DB model with Te//=18621 K, log(g) 
= 8.2, d = 85 pc, Rwd = 0.009 R©. The resultant model is plotted as a dashed grey line, 
while the WD model with the background galaxy added is shown as a dash-dot-dot-dot line. 
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The models for SDSS0845 are shown in Figure[2l The first shows the best-fitting disk 
model, with an inner temperature of 1800 K, an outer temperature of 250 K, and an inclina- 
tion of 85 degrees. Figurel^b shows the closest physical model with an inner disk temperature 
constrained to 1400 K, Tout = 400 K and inclination = 80 degrees. This corresponds with 
inner and outer disk radii of Rm = 18.8 Rwd {0.17 Rq) and Rout = 99.7 Rwd (0.89 7?©). 
The coolest disk model still within 3a of the plotted datapoints has = HOOK, Tout = 
600 K and an inclination of 70 degrees. 

The outer radius of the gaseous component was estimated by Gaensicke et al. (2008) to 
be ~ 0.8 — 0.9Rq, which is comparable to that of the dust disk derived above. 



3.3. SDSS J1043+0855 



The GALEX, UKIDSS and Spitzer flux densities for SDSS1043 are given in Table^l 
and are plotted in Figure[3] together with the SDSS DR7 spectrum. The uncertainties on the 
Spitzer photometry for SDSS1043 are notably large, due to the faint nature of the source 
(only 30 fiJj in IRAC channel 1). In addit i on to these data points, we include CFHT and 
Lick near-IR photometry from iMelis et al.l (120101 ) to constrain our model. All previously- 
published data are shown in dark grey in Figure[3l while the GALEX, UKIDSS and Spitzer 
data are shown in light grey. We hav e fitted the DR 7 spectrum of SDSS1043 using DA model 
spectra computed with th e code of Koesteil (2010), w hich incorporates the updated Stark 
broadening calculations of iTremblay fc Bergeron! (120091 ) . and find Tg// = 17912 ± 362 K and 
log{g) = 8.07±0.08, corresponding to Mh^z) = O.66±O.O5M0 and Rwd = 8.67±0.47x 10^ cm 
(Table 3). These v alues are, within the un certainties, consistent with those derived from the 
DR4 spectrum by iGansicke et al.l (120071 ). but given the improvement in the SDSS data 
reduction and line broadening, the new values reported here should be preferred. 

The IRAC data points show a 3cr excess above the expected flux density of the white 
dwarf, increasing dramatically at 8 microns. 

As with the other targets, we show the best-fit optically thick dust disk model in Fig- 
ure[3^, with an inner disk dust temperature of Tj„ = 1800K, Tout = 1700 K and an inclination 
of 40 degrees. When constraining the inner disk temperature below the sublimation thresh- 
old of 1400 K, we find a best-fit model with Tj„ = 1400, Tout = 800 K and i = 85 degrees. 
This corresponds with inner and outer disk radii of Rin = 17.9 Rwd (O.23-R0)and Rout = 
37.6 Rwd (0.497?©). The coolest model allowed by the fit has T^^ = 900 K, Tout = 800 K and 
i = 56 degrees. 



Modelling the double-peaked Call line profiles. 



Gansicke et al, 



( 120121 ) find, similar as 
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for SDSS0738 and SDSS0845, 



IRq (see also lMelis et al.ll2010h 



3.4. SDSS J1617+1620 

The flux densities from GALEX, SOFI and Spitzer are given in TableEj and are plotted 
in light grey in FigurelH together with the SDSS DR7 spectrum (dark grey). The data show 
a clear excess above the flux den sity from the white dw arf at the Spitzer wavelengths. The 
white dwarf model, also based on lGansicke et al.l ( 120111 ). is shown as a dashed black line and 
is a DA WD model, with T^// = 13432 K, log(g) = 8.04, d = 122 pc, Rwd = 0.01261 R©, 
reddened by E(B-V) = 0.04 to better match the optical SDSS spectrum. The reddening has 
negligible effect at infrared wavelengths. 

The best-fit dust disk model has an inner disk effective temperature Tj„ = 1800 K, 
outer disk Tout = 950 K and an inclination of 70 degrees, plotted as a grey dash-dot line 
in Figure!!^. The closest fit with Tj„ fixed at 1400K has Tout = 950 K and i = 50 degrees 
(FigureHb). This corresponds with inner and outer disk radii of Rm = 12.2 Rwd {0.15Rq) 
and Rout = 20.4 Rwd (0.26 i?©). The model with the coolest possible inner disk temperature 
has Tin = HOOK, Tout = 620 K and i = 53 degrees. 

The peak separation of the Call emission lines i mplies a maximum outer radius of the 
gaseous component of 90Rwd (IGansicke et al.l 120121 ). 



Discussion and Conclusions 



We find that out of these four white dwarfs with previously- identified gaseous disks, 
all show evidence for dusty extensions to those disks, as seen in the previously-studied 
SDSS J1228+1040. We note that the best-fitting models for all of these disks have inner disk 
temperatures hotter than the sublimation temperature expected from the bulk composition 
of the dust, although all can be modelled satisfactorily with inner disk temperatures fixed 
to a physically more plausibl e temperature of 1400 K . When adding in the results from our 
previous study of SDSS1228 (IBrinkworth et al.ll2009l ). we find that all five white dwarfs with 
Call emission from a gaseous debris disk also show dusty components to those disks. It 
is evident that gaseous and dusty debris disks easily c oexist aroun d wh ite dwarfs, a r id the 
dusty disks likely feed the gaseous component. Both lJural (120081 ) and iMelis et al.l (120101 ) 
suggest that dust disks are a pre-requisite to gaseous disks, and our findings are consistent 
with that evolutionary path, and vice- versa, the presence of gas is likel y to have a ma jor 
effect on the evolution of the planetary debris discs around white dwarfs (jRafikovll2011bl ). 
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The one outlier of our sample is SDSS1043. Even taking the uncertainties on the 
photometry into account, it is apparent that the dust in the system is far less pronounced 
than in the other three targets, despite very similar gaseous disk signatures in all four 
systems. The 4.5/im excesses seen in SDSS0738, SDSS0845 and SDSS1617 are 7.5, 5.7, and 
4.5 times the white dwarf contribution, respectively. In contrast, in SDSS1043, the measured 
flux density at 4.5 microns is only 1.7 times the white dwarf contribution. Why SDSS1043 
should lack significant dust compared to the other three systems remains a mystery. Maybe 
the system is at a higher inclination, so we are seeing the correspondingly smaller projected 
area of the edge of the geometrically thin disk, rather than its face, but this doesn't seem 
to be borne out in the modelling of the other systems. Considering the likelihood that 
the dust disks feed the gaseous disks in these systems, it is feasible that we are seeing the 
tail-end of an asteroidal tidal disruption event that took place at an earlier time than the 
other systems, such that the dust around SDSS1043 is now depleted compared to the other 
three white dwarfs, and the gaseous disk dominates the SED. This seems consistent with the 
most recent estimates for disk lifetimes: a massive (~ 10^^) disk accreting due to Poynting- 
Robertson drag alone can be expected to survive for several Myr ( iRaf okovl |2 1 1 al ) . However, 
when a gas disk is present i r i the s ystem due to sublimation of the inner disk, runaway 
accretion can occur. iRafikov J2011bh predict that this can lead to the disk being exhausted 
in ~ 10^ yr. iBochkarev &: Rafikovl (l201l[ ) extend this analysis to include gas generated from 
multiple asteroid in ipacts, i n whic h they confirm the estimated disk lifetime of several Myr. 
On the other hand, iDebesI ( 1201 Ibl ) calculates that dust disks around white dwarfs may last 
for up to 1 Gyr, while gas disks that are not being replenished are unlikely to last more 
than 10^ years before fully accreting onto the white dwarf. If the disk around SDSS1043 is 
indeed reaching the end of its lifetime, long-term infrared monitoring of the SEDs of all of 
the five systems (SDSS1228, SDSS0738, SDSS0845, SDSS1043 and SDSS1617) over the next 
few decades may help to shed some light on this. 



The authors would like to thank the referee for detailed and useful comments that helped 
to improve this paper. 

The DA white dw arf models used in the analysi s were calculated with the modified Stark 



broadening profiles of iTremblay fc Bergeron! ( l2009l ). kindly made available by the authors. 



Co-authors Gaensicke and Marsh were supported by an STFC rolling grant. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, Caltech, under NASA contracts 1407 and 960785. 
This work makes use of data products from the Two Micron All Sky Survey, which is a joint 
project of the University of Massachusetts and IPAC/Caltech, funded by NASA and the 



- 11 - 



NSF. The SOFI results in this paper are based on observations collected at the European 
Southern Observatory (La Silla) under programme ID 085.D-0541. This publication makes 
use of data products from the Wide-field Infrared Survey Explorer, which is a joint project 
of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California 
Institute of Technology, funded by the National Aeronautics and Space Administration. 
Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the 
Participating Institutions, the National Science Foundation, the U. S. Department of Energy, 
the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the 
Max Planck Society and the Higher Education Funding Council for England. The SDSS Web 
Site is |http: / /www.sdss.or^/] The SDSS is managed by a very long list of institutions, who 
can all be found at the following website: http:/ /www. sdss.org/collaboration/credits. html. 

This research has made use of NASA's Astrophysics Data System. This research has 
made use of the SIMBAD database, operated at CDS, Strasbourg, France. 

REFERENCES 

Abazajian et al. 2009, ApJ, 182, 543 

Adelman-McCarthy, J.K., Agiieros, M.A., Allam, S. S., Anderson, K. S.J., Anderson, S.F., 
Annis, J., Bahcall, N. A., Baldry, I. K., Barentine, J. C, Berlind, A. and 131 coauthors, 
2006, ApJS, 162, 38 

Becklin, E. E., Farihi, J., Jura, M., Song, I., Weinberger, A. J., &; Zuckerman, B. 2005, ApJ, 
632, L119 

Bertin, E., Arnouts, S., 1996, A&AS, 117, 393 
Bochkarev, K. V., Rafikov, R. 2011, ApJ, 741, 36 

Brinkworth, C. S., Gansicke, B. T., Marsh, T. R., Hoard, D. W., Tappert, C. 2009, ApJ, 
696, 1402 

Brinkworth, C. S., Hoard, D.W., Wachter, S., Howell, S.B., Ciardi, D.R., Szkody, P., Har- 
rison, T.E., van Belle, G.T., Esin, A. A. 2007, ApJ, 659, 1541 

Debes, J. H., Hoard, D. W., Kilic, M., Wachter, S., Leisawitz, D. T., Cohen, M., Kirkpatrick, 
J. D., Griffith, R. L. 2011, ApJ, 729, 4 

Debes, J. 2011b, White Dwarf Atmospheres and Circumstellar Environments (Ed. D. W. 
Hoard), Wiley. 



- 12 - 



Debes, J.H., Sigurdsson, S. 2002, ApJ, 572, 556 

Dufour, P., Kilic, M., Fontaine, G., Bergeron, P., Lachapelle, F, -R., Kleinman, S. J., Leggett, 
S. K., 2010, ApJ, 719, 803 

Dufour, P., Kilic, M., Fontaine, G., Bergeron, P., Melis, C., Bochanski, J., 2012, ApJ, 
accepted. larXiv:1201.6252i 

Fazio, G., et al. 2005, ApJS, 154, 10 

Farihi, J., Becklin, E. E., Zuckerman, B., 2005, ApJS, 161, 394 
Farihi, J., Zuckerman, B., Becklin, E.E., 2008, ApJ, 674, 431 
Farihi, J., Jura, M., Zuckerman, B. 2009, ApJ, 694, 805 
Farihi, J., Jura, M., Lee, J.-E., Zuckerman, B. 2010, ApJ, 714, 1386 

Ferguson, H. C., Dickinson, M., Giavalisco, M., Kretchmer, C., Raavindranath, S., Idzi, 
R., Taypor, E., Conselice, C.J., Fall, S.M., Gardner, J. P., Livion, M., Madau, P., 
Moustakas, L, A., Papovich, C. M., Somerville, R. S., Spinrad, H., Stern, D. 2004, 
ApJL, 600, L107 

Gansicke, B. T., Marsh, T.R., Southworth, J., Rebassa-Mansergas, A. 2006, Science, 314, 
1908 

Gansicke, B.T., Marsh, T.R., Southworth, J. 2007, MNRAS, 380, L35 

Gansicke, B. T., Koester, D., Marsh. T. R., Rebassa-Mansergas, A., Southworth, J. 2008, 
MNRAS, 391L, 103 

Gansicke, B. T. 2011, in Planetary systems beyond the main sequence, AlP Conf. Proc. 1331, 
211 (Eds. S. Schuh, H. Drechsel,& U. Heber) 

Gansicke, B. T., et al. (2012) in prep. 

Girven, J., Brinkworth, C. S., Farihi, J., Gansicke, B. T., Hoard, D. W., Marsh, T. R., 
Koester, D. 2012, accepted ApJ. larXiv:1202.3784V l 

Graham, J.R., Matthews, K., Neugebauer, G., Soifer, B. T., 1990, ApJ, 357, 216 

Hambly, N. C., Collins, R.S., Cross, N.J.G., Mann, R.G., Read, M.A., Sutorius, E.T.W., 
Bond, I., Bryant, J., Emerson, J. P., Lawrence, A. and 7 co-authors. 008, MNRAS, 
384, 637 



- 13 - 



Hewett, R. J., Warren, S.J., Leggett, S.K., Hodgkin, S.T. 2006, MNRAS, 367, 454 

Hoard, D.W., Howell, S.B., Brinkworth, C.S., Ciardi, D.R., Wachter, S. 2007, ApJ, 671, 
734 

Houck, J. R., et al. 2004, ApJS, 154, 211 
Jura, M. 2003, ApJ, 584, L91 

Jura, M., Farihi, J., Zuckerman, B. 2007, ApJ, 663, 1285 
Jura, M., AJ, 135, 1785 

Kilic, M., von Hippel, T., Leggett, S.K., Winget, D. E. 2005, ApJ, 632, L115 
Kilic, M., von Hippel, T., Leggett, S.K., Winget, D. E. 2006, ApJ, 646, 474 
Kilic, M., Redfield, S. 2006, ApJ, 660, 641 

Kilic, M., Patterson, A. J., Barber, S., Leggett, S. K., Dufour, R 2011, MNRAS, 419, L59 
Klein, B., Jura, M., Koester, D., Zuckerman, B. 2011, ApJ, 741, 64K 
Koester, D., Provencal, J., Shipman, H. L. 1997, A&A, 320, L57 
Koester, D. 2008, arXiv0812.0482 
Koester, D. 2010, MmSAl, 81, 921 

Kuchner, M.J., Koresko, CD., Brown, M.E., 1998, ApJ, 508, L81 

Lawrence, A., Warren, S. J., Almaini, O., Edge, A.C., Hambly, N.C, Jameson, R. P., Lucas, 
P., Casili, M., Adamson, A., Dye, S. and 12 co-authors. 2007, MNRAS, 379, 1599 

Lodders. K. 2003, ApJ, 591, 1220 

Makovoz, D., Roby, T., Khan, I., Booth, H., 2006, SPIE, 6274, 10 
Maraston, C. 2005, MNRAS, 362, 799 
Maraston, C. 1998, MNRAS, 300, 872 
Martin et al. 2005, ApJ, 619, 1 

Melis, C, Jura, M., Albert, L., Klein, B., Zuckerman, B. 2010, ApJ, 722, 1078 
Patten, B., et al. 2006, ApJ, 651, 502 



- 14 - 



Rafikov, R. 2011, ApJL, 732, L3 
Rafikov, R. 2011, MNRAS, 416L, 55 

Reach, W. T., Kuchner, M. J., von Hippel, T., Burrows, A., Mullally, F., Kilic, M., Winget, 
D.E. 2005, ApJ, 635, L161 

Rieke, G. H. et al. 2004, ApJS, 154, 25 

Rousselot, R, Lidman, C, Cuby, J.-G., Moreels, G., Monnet, G. 2000, A&A, 354, 1134 

Skrutskie et al. 2006, A J, 131, 1163 

Telesco, C. M., Joy, M., Sisk, C., 1990, ApJL, 358, 17 

Tokunaga, A. T., Becklin, E. E., Zuckerman, B., 1990, ApJL, 358, 21 

Tremblay, R-E., Bergeron, R 2009, ApJ, 696,1755 

Vennes, S. Kawka, A., Nemeth, R 2011, MNRAS, 413, 2545 

von Hippel, T., Kuchner, M. J., Kihc, M., Mullally, P., Reach, W. T. 2007, ApJ, 662, 544 
Xu, S., Jura, M., 2011, ApJ, in press. l arXiv:1109.4"207l 
Zuckerman, B., & Becklin, E. E. 1987, Nat., 330, 138 

Zuckerman, B., Koester, D., Mehs, C., Hansen, B. M., Jura, M. 2007, ApJ, 671, 872 



This preprint was prepared with the AAS E^TJrjX macros v5.2. 



- 15 - 



Table 1. Observation log for Spitzer Space Telescope data 



Target 


Instrument/ 


AOR Key 


Total Integration 


Date 


Pipeline 




Channel 




Time (s) 






SDSS J0845+2257 


IRAC Ch. 2 


25459968 


600 


2008 May 09 


S18.7 


SDSS J0845+2257 


IRAC Ch. 4 


25459968 


600 


2008 May 09 


S18.7 


SDSS J0845+2257 


IRS Blue PUI 


25460224 


600 


2008 Dec 02 


S18.7 


SDSS J0845+2257 


MIPS- 24 


25460480 


140 


2008 May 18 


S18.12 


SDSS J0845+2257 


MIPS-24 


25460736 


140 


2008 May 19 


S18.12 


SDSS J1043+0855 


IRAC Ch. 1 


22248448 


3000 


2008 Jun 10 


S18.7 


SDSS J1043+0855 


IRAC Ch. 2 


22248448 


3000 


2008 Jun 10 


S18.7 


SDSS J1043+0855 


IRAC Ch. 3 


22248448 


3000 


2008 Jun 10 


S18.7 


SDSS J1043+0855 


IRAC Ch. 4 


22248448 


3000 


2008 Jun 10 


S18.7 


SDSS J1617+1620 


IRAC Ch. 1 


39873024 


750 


2010 Sep 01 


S18.18 


SDSS J1617+1620 


IRAC Ch. 2 


39873024 


750 


2010 Sep 01 


S18.18 


SDSS J0738+1835 


IRAC Ch. 1 


39872768 


3000 


2010 Dec 01 


S18.18 


SDSS J0738+1835 


IRAC Ch. 2 


39872768 


3000 


2010 Dec 01 


S18.18 



Table 2. Photometry for all four target stars. 



V V dVCUclllLl 


lA/QTT'olon (T'J' n { //TYl 1 

VVclVddlg,Lll I ^lli J 




Flux Density (mJy) 








SDSS0738+1835 


SDSS0845+2257 


SDSS1 043+0855 


SDSS1617+1620 

k_/ J— ' k_/ k_/ J- \_/ J- 1 1 £^\J 


GALEX FUV 


0.158 


0.061+}^P 


1.92±5% 


0.66±5% 


0.345+^4! 

— 4.4/0 


GALEX NUV 


0.227 


0.161+tI? 

—4.970 


2.28±5% 


0.54±5% 


0.406+^3^^° 

— 1.070 


Y 


1.02 







0.182+^-9%(6) 




J 


1.2 


_ 


0.473±ll%(a) 


0.128+^-^^°(6) 


_ 


H 


1.6 


_ 


0.349±14%(a) 


0.088l^:g(6) 




K 


2.2 


- 


0.294±17%(a) 


0.049+S(6) 


0.122^-5(c) 


WISE-1 


3.35 




0.267±4.5% 




0.104±7.7% 


IRAC-1 


3.550 


0.082±3% 




0.039±6% 


0.108±3% 


IRAC-2 


4.493 


0.075±3.7% 


0.248+^% 


0.025±8% 


0.095±3% 


WISE-2 


4.60 




0.237±8.9% 




0.097±15.5% 


IRAC-3 


5.731 






0.017±16% 




IRAC-4 


7.872 




n 1 70+21% 


0.029±27% 




IRS PU-16 


15.60 




n ic; 1+3.6% 
^■-'■^■'■-4.7% 






MIPS-24 


23.68 




U--^J--^_4% 







Note. — Source of NIR photometry: (a) AAT (b) UKIDSS (c) SOFT 
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Table 3. Stellar parameters for target stars 



Parameter 


SDSS J0738+1835 


SDSS J0845+2257 


SDSS J1043+0855 


SDSS J1617+1620 


WD Temperature (K) 


13600 


18600 


17912 


13432 


Distance (pc) 


120 


85 


183 


122 


Iog{g) 


8.5 


8.2 


8.0 


8.0 


Radius (cm) 


6.16e8 


6.27e8 


9.05e8 


8.77e8 





Wavelength (microns) 



Wavelength (microns) 



Fig. 1. — The spectral energy distribution (SED) of SDSS J0738+1835 from the ultraviolet 
to the mid-infrared. See Table|2]for the GALEX and Spitzer data. The WD spectrum was 
taken from the SDSS Data Release 7 and is plotte d along with t h e GA LEX, Spitzer, WISE 
and SDSS photometry, and the near-IR data from iDufour et al.l (|2010| ). The SED shows a 
significant flux density excess above that expected from the white dwarf alone. The SED 
on the left is compared to a white dwarf + optically thick disk model with an inner dust 
temperature of 1800 K, and outer disk temperature of 830 K and an inclination of 58 degrees. 
The SED on the right limits the inner disk temperature to Tj„ = 1400 K, with Tout = 930 K 
and a face-on inclination of degrees. See Section [3?T] for the full range of possible models. 
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Fig. 2.— The SED for SDSS J0845+2257 from the ultraviolet to the mid-infrared GALEX, 
AAT and Spitzer d ata are plotted in light grey, while SDSS, NIR (iMelis et al.lboiol ). AKARI 
(iFarihi et al.ll2010l ) and WISE Preliminary Data Release photometry are plotted in dark grey, 
along with the SDSS DR? white dwa r f spec t rum. The white dw arf model plus background 
galaxy (8 Gyr, z=1.6. iMaraston et al.l ( 119981 ) : lMaraston. C.I (120051 ): see Section [273|) is plotted 
in dash-dot-dot-dot. The SED shows a significant flux density excess from the K-band to 
longer wavelengths, over the flux density expected from the white dwarf and the background 
galaxy alone. The SED on the left is compared with the best model of a WD -|- background 
galaxy -|- optically thick dust disk with an inner dust temperature of 1800 K and an outer 
disk temperature of 250 K. The disk inclination is 83 degrees. We note that, while this is a 
good fit, the inner dust temperature is above the expected dust sublimation temperature. 
The SED on the right is compared to the same model but with a physically-possible inner 
dust disk temperature of 1400 K, an outer disk temperature of 400 K and an inclination of 
80 degrees. See Sectioning] for the full range of possible model parameters. 



- 19 - 




0.001 I 

0-1 



Wavelength (microns) 




Wavelength (microns) 



Fig. 3.— The SED of SDSS 1043+0855 from the ultraviolet to the mid- infrared. The WD 
spectrum was taken from SDSS Data Release 7 and is plotted along with t he GALEX 



UKID SS and Spitzer in light grey, and the SDSS photometry and NIR data from lMelis et al. 



(I2OIOI ) in dark grey. The SED shows a significant flux density excess above that expected from 
the white dwarf alone. The SED on the left is compared to a white dwarf + optically thick 
disk model with an inner dust temperature of 1800 K, and outer disk temperature of 1700 K 
and an inclination of 40 degrees. The SED on the right limits the inner disk temperature to 
Tin = 1400 K, with Tout = 800 K and an inclination of 85 degrees. See Section [373] for the full 
range of possible models. 
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Fig. 4. — The SED for SDSS J1617+1620 from the ultraviolet to the mid-infrared. Spitzer, 
GALEX, WISE and SOFT data are reported in TableU GALEX, SOFT and Spitzer data are 
plotted in light grey, while the SDSS data and the WISE Preliminary Data Release are shown 
in dark grey. The WD spectrum was taken from SDSS Data Release 7. The SED shows a 
significant flux density excess from the K-band to longer wavelengths, over the flux density 
expected from the white dwarf alone. The SED on the left is compared to the best-fit white 
dwarf -|- optically thick disk model with an inner dust temperature of 1800 K, and outer disk 
temperature of 950 K and an inclination of 70 degrees. The SED on the right limits Tj„ to 
1400 K with Tout = 950 K and i = 50 degrees. See Section l3^ for the full range of model 
parameters. 



